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Abstract 



We address the problem of constructing positive operator-valued 
measures (POVMs) in finite dimension n consisting of n 2 operators 
of rank one which have an inner product close to uniform. This is 
motivated by the related question of constructing symmetric informa- 
tionally complete POVMs (SIC-POVMs) for which the inner products 
are perfectly uniform. However, SIC-POVMs are notoriously hard to 
construct and despite some success of constructing them numerically, 
there is no analytic construction known. We present two constructions 
of approximate versions of SIC-POVMs, where a small deviation from 
uniformity of the inner products is allowed. The first construction 
is based on selecting vectors from a maximal collection of mutually 
unbiased bases and works whenever the dimension of the system is 
a prime power. The second construction is based on perturbing the 
matrix elements of a subset of mutually unbiased bases. 

Moreover, we construct vector systems in C n which are almost 
orthogonal and which might turn out to be useful for quantum com- 
putation. Our constructions are based on results of analytic number 
theory. 

1 Introduction 

1.1 Background 

A bàsic question in quantum mechanics is how to obtain information about 
the state of a given physical system by using suitable measurements. Even 
in case many identically prepared copies of the system are available, it is a 
nontrivial task to devise a measurement procedure which uniquely identifies 
the given quantum state from the statistical data produced by the measure- 
ments. Note that this holds true even in case the complete statistics, that is, 
the probabilities for the different measurement outcomes, is known. 

We next describe the possible measurements of the quantum system in 
more detail and first remark that all systems considered in this paper are of 
finite dimension n. If the state of the quantum system is given by an n x n 
density matrix, then the complete measurement statistics of one fixed von 
Neumann measurement is not suficient to reconstruct the state. Indeed, this 
follows from the fact that the statistics of a fixed von Neumann measurement 
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determines at most n — 1 real parameters (specified by the probabilities of 
the measurement outcomes), whereas a general density matrix is determined 
by n 2 — 1 free real parameters. 

In fact it is possible to perform a more general measurement procedure 
on a quantum system, namely a posítive operator-valued measure, or POVM 
for short, see j33J HH]. A POVM is described by a collection of positive 
operators Ei > 0, called POVM elements, that partition the identity, that 
is, Y^í-Eí = I n - If the state of the quantum system is given by the density 
matrix p, then the probability Pi to observe outcome i in the POVM is given 
by the Born rule 



where ti{A) denotes the trace of a complex matrix A. The task is to devise 
a POVM with operators Ei such that the state p is uniquely specified by 
the probabilities Pi in (0). The POVM is then called ínformatíonally com- 
plete, or simply an IC-POVM, and they appear to have been first studied 
in [SE]. A particularly interesting question is whether a POVM exists on 
C n that consists of n 2 POVM elements Ei of rank one. Counting the num- 
ber of parameters determined by the measurement, we see that n 2 is indeed 
the minimal possible number of such POVM elements. In this case Ei is 
a subnormalized projector, that is, Ei = Hi/n for projectors Ilj = \ipi) (^| 
corresponding to some vectors \ipi) in C n . In |T2] it has been shown that 
IC-POVMs exist in all dimensions and in ^3] a method has been given how 
to construct IC-POVMs by taking a fixed fiducial start vector and taking 
the orbit of this vector under a (projective) group operation. 

As an example of this type, consider the normalized states j^i) , • • ■ , \ip4) 
defined as follows: 



Then the rank one operators defined by Ei = 1/2 1^) (ijji\ are given by 



Pi = ti(pEi) 



(1) 
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and it can be verified easily that E\ + E 2 + E3 + E4 = J4 is the identity 
matrix and that the matrices E\, E 2 , E 3 , E4 are linear ly independent. For 
the possible inner products between two POVM elements Ei and Ej where 
i^jwe obtain that ti(EiEj) G {4/81,49/324}. 

Our goal is to find IC-POVMs on C n such that nhi(EiEj) is "small" for 
distinct POVM elements Ei = \ipi)(ipi\/n and Ej = \ipj)(ipj\/n. In Section fOl 
we make precise what we mean by the inner products being small and briefly 
summarize previous work on the problem. 

1.2 Notation 

We use the Landau notation to compare the asymptotics of two functions 
f,g : N — > C. We recali that /(n) = o(g(n)) means lim^oo f(n)/g(n) = 0. 
Furthermore, f(n) = 0(g(n)) means that there exists a constant c > 0, 
such that \f(n)\ < cg(n) for all n > 1. Throughout the paper, the implied 
constants in the symbols 'o' and 'O' may occasionally, where obvious, depend 
on an integer parameter d and a small real parameter e > 0, and are absolute 
otherwise. 

For an integer n > 1 we denote by C n the n-dimensional vector space 
over the complex numbers C. For two vectors = (a\, . . . , a n ) G C n and 
\<p) = (61, . . . , b n ) G C n , we use 

n 
i=l 

to denote their inner product. We also define 5 it j — 1 if i — j and otherwise. 
We denote the identity matrix of size n x n by I n , the all-ones matrix of size 
n by J n = [l]ij =1 - We use diag(ai, . . . ,a n ) to denote the n x n diagonal 
matrix which has a 1; . . . , a n on the main diagonal. For matrices A and B, 
we use A © B to denote their block-diagonal direct sum. 

For a real z and an integer m we use the notation e m (z) = exp(2iriz/m), 
where 1 = 

We use Fq to denote the finite field of q elements, and we also assume 
that for a prime p, the field ¥ p is represented by the set {0, . . . ,p — 1}. 

As we have mentioned, we use tr(v4) to denote the trace of a complex 
matrix A. On the other hand, for an element a G ¥ q we use Trf q /f p (a) to 
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denote its trace in the prime subfield ¥ p of F 9 , see jSI]. That is, if q = p m 
for a prime m, then 

m—X 

Tr F 9 /F p (a) = 5^ aPJ - 

3=0 

1.3 Previously Known Results 

A particularly appealing case of IC-POVMs arises when furthermore all of 
the inner products of the vectors are small. An extremal case in this 
sense arises when we are given a system of n 2 normalized vectors {\if>i) '■ i = 
1 , . . . , n 2 } in C n for which 

m\^)? = ^l^ l<K]<n 2 . (2) 

Indeed, for any system of k vectors j^i) , • • • , \ipk) i n C n for which the absolute 
vàlues of pairwise inner products are constant | (V^l^j) | 2 = oc, where oc G R 
(for 1 < i < j < k) the so-called special bound holds which says that 
k < "fi ■ Specializing a = l/(n + 1) we obtain that n 2 is the largest 
possible number of vectors satisfying (J2J). 

A system of vectors as in (J2J), respectively the corresponding POVMs, 
are called symmetríc informationally complete POVMs, or SIC-POVMs for 
short. They have se ver al very desirable prop ert ies, see ^2] for a discussion 
in the context of the quantum de Finetti theorem and more generally in 
their Bayesian approach to quantum mechanics and its interpretation [TT] . 
Furthermore, see |THJ QÜ] for their role in establishing the quantumness of 
a Hilbert space and the related question about optimal intercept-resend 
eavesdropping attacks on quantum cryptographic schemes. Explicit ana- 
lytical constructions of sets satisfying (J2J have been given for dimensions 
n = 2,3,4,5 in [3Ü|. Section 3.4] and [41) . for dimension n = 6 in [20], for 
dimensions n = 7, 19 in £Q and for dimension n = 8 see [21]. While it has 
been conjectured that SIC-POVMs exist in all dimensions jSÜJ Section 3.4] 
or jJT] and numerical evidence exists for dimensions up to 45, see it is 
a difficult task to explicitly construct systems of vectors which satisfy (j2J). 
There are no known infmite families of SIC-POVMs and in fact it is not even 
clear if there are SIC-POVMs for infmitely many n. 
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1.4 Our Results 



In the first part (Section |2J of this paper we relax condition (J2J) on the inner 
products slightly and allow that 

l(^i)| a <— , \<i<j<n 2 . (3) 

The purpose of the first part of this paper is to show that infinite fam- 
ílies of systems of n 2 normalized vectors which satisfy Q and give rise to 
IC-POVMs exist. We call the rank one projectors obtained from such sys- 
tems of vectors approximately symmetríc informationally complete POVMs, 
or ASIC-POVMs for short. Here we show that when n = p r is a power of a 
prime p, ASIC-POVMs can be constructed. 

In the second part (Section 0} of the paper we explore properties of other 
approximately symmetric vector systems where we do not require the proper- 
ties of completeness and informational completeness but require the property 
of approximate symmetry. We also relax (JHJ) further by allowing that the in- 
ner products be bounded from above by \{ipí\ipj)\ < (2 + o(l))/y/n and by 
dropping the requirement that the vectors give rise to a POVM. This addi- 
tional freedom then allows us to construct bases in all dimensions n. Besides 
their general mathematical interest in constructing such vector systems they 
might be useful in quantum cryptographic scenarios which generalize the 
BB84 setting |7j, such as the protocols described in |H1 E2H Eiïj • See also [TJ] 
for an analysis of general schemes for quantum key distribution where the 
sender uses arbitrary quantum states and the receiver's measurement is re- 
placed by a POVM. 

Besides approximations to SIC-POVMs we also consider approximations 
to mutually unbiased bases (MUBs). Since we also need MUBs for our con- 
struction of ASIC-POVMs we briefly recali their definition. A maximal set 
of MUBs is given by a set of n 2 + n vectors in C n which are the elements of 
n + 1 orthonormal bases = {l^fc,].) , • • • , \ipk,n)} °f C n where k — 0, . . . , n. 
Hence, 

(V'JM^fcj) = ( 4 ) 
and the defining property is the mutual unbiasedness, given by 

KVjmIV*j>I = 4= (5) 
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for < k,£ < n, k ^ £, and 1 < i,j < n. Starting from fM EE21 liü] an 
extensive growing body of research explores MUBs and their constructions, 
see E El 113 EE3 Hi] 123 ES EU EH] and references therein. However, so 
far maximally sets of n + 1 MUBs in dimension n are only known to exist 
in any dimension n = p r which is a power of a prime p > 3, see [27] for an 
overview and some of such constructions. The main construction is based on 
Gaussian sums and in the case of prime n = p can be described as 

\ip k ,j) = (e p (ku 2 + ju)) p u=1 , 1 < k,j <p, 

and also £> being a Standard orthonormal basis, that is, \ipoj) = (<^>)u =r 

One can use additive characters over an arbitrary finite field to extend 
this construction to an arbitrary prime power n = p r . However the condition 
that n = p r is a prime power is still somewhat too restrictive and unnatural 
for quantum computation. So a natural question arises whether MUBs exist 
for every positive integer n. In the second part we consider vector systems 
where we relax the conditions © and (JHJ). We use exponential sums to 
construct vector systems for any dimension n which 

• satisfy (jlj) but instead of © all other inner products are 0(n~ 1//4 ); 

• is normalized but instead of (0) all other inner products are at most 
(2 + o(l))n~ 1 / 2 . 

We call vector systems of n 2 + n vectors in C n which satisfy (@J), and 
instead of (jïïj) the condition 

< 1<K] <n 2 + n 

approximately mutually unbiased bases, or AMUBs for short. 

We also construct some vector systems using multiplicative and mixed 
character sums, which 

• satisfy (@J) and assuming some natural and widely believed conjecture 
on the distribution of primes in arithmetic progressions all other inner 
products are 0(n _1//2 log?7,); 
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• in the special case of n = p — 1 where p is a prime, form AMUBs. 

Interestingly, our arguments use both the classical bound of Weyl |27] 
(see also [23 E3) as well as the more recent, but no less celebrated, bounds 
of Weil (see also JHül EU [46 ) . Besides exponential sum techniques we 
also use some recent results about the gaps between prime numbers from [3] . 
We conclude with some conjectures and open qüestions concerning our con- 
structions in SectionHJ 

2 Constructing ASIC-POVMs 
2.1 Preliminaries 

We begin by giving a defmition of the vectors and associated rank one oper- 
ators we are interested in. 

Definition 1 (ASIC-POVMs). Suppose that n is a posítive integer. Let 
•A = (IV'í) : i = l,---,"^ 2 } be a set of vectors in C n . Let £ = {Ei = 
\ipi) /n '■ i — 1, ■ ■ ■ , n 2 } be the corresponding set of subnormalized pro- 
jection operators. If S satisfies the condítions 

2 

(i) Ym=\Eí = I n fcompleteness/POVM conditionj; 

(zz) the matrices Ei are linearly independent as elements ofC nxn (Informa- 
tional completenessj; 

{in) \{ipi\ipj)\ 2 = nHr(EiEj) < (1 + o(l))n _1 for 1 < i < j < n 2 (approxi- 
mate symmetryj; 

then we call S an approximately symmetric informationally complete POVM, 
or ASIC-POVM for short. 

We remark, that in fact, sometimes we also refer to the corresponding set 
A as an ASIC-POVM. 

In the subsequent sections, we present two different constructions that 
give rise to infinite families of ASIC-POVMs. The first construction is based 
on the observation that a set of n + 1 mutually unbiased bases in C n gives 
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rise to an IC-POVM, cf. j^ü- However, this IC-POVM consists of n 2 + n 
rank-one operators; thus, it is an overcomplete generating set of the vector 
space of all n x n matrices. In our first construction in Section 12.31 we show 
how to select n 2 projectors that allow us to derive an ASIC-POVM. The 
second construction in Section 12.41 starts from all vectors contained in n of 
the n + 1 MUBs. We show that by slightly perturbing the components of 
these vectors it is possible to obtain an ASIC-POVM. 

2.2 POVMs and Frames 

Suppose that A = G C n : 1 < i < n 2 } is a system of n 2 vectors of unit 

norm, such that A spans C n and the associated subnormalized projectors 
Ei = \ /n satisfy n 2 tr(EiEj) = (1 + o(l))ri _1 whenever i ^ j. 

We would like to have that the subnormalized projectors Ei form a 

2 

POVM, but, unfortunately, the completeness relation Ym=i Ei = ls m 
general not satisfied. However, there is a way to fix this using a technique 
described in : Define a positive semidefinite hermitian operator G by 

n 2 

1=1 

Since A spans C ra , the inequality (<p\G\<p) = J2i li^il^l 2 > holds for any 
nonzero vector |<£>) in C n , so G is even positive definite. It follows that G^ 1 
exists and is positive definite, and we can form the uniquely determined 
positive definite square-root G~ x l 2 . The n 2 rank-one operators 

£ — {Fi — G-^E.G- 1 ' 2 : \<i<n 2 } 

form a POVM, since £\ F t = G^^GG^I 2 = I n . 

Clearly, if the operators Ei are linearly independent, then so are the 
operators Fi. Therefore, the procedure preserves information-completeness. 

In general, if we switch from the rank-one operators Ei to the rank-one 
operators Fi, then tr(FiFj) < (1 + o(l))/n 3 might not hold for some i ^ j. 
However, if G^ 1 is close to the identity matrix, then approximate symmetry 
is preserved as well. 

We now mention some connections between POVMs and the theory of 
frames 00 US]. 
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Definition 2 (Frames). A set T = {\ipi) : 1 < i < N} of vectors in C n is 
called a frame if there exist real numbers a and b, with < a < b, such that 

N 
i=l 

holdsfor all \<p) E C n . 

If a = b, then the frame is called a tight frame, and if a = b = 1 then the 
frame is called a Parseval frame. 

We can associate with the frame T its frame operator G = Ylk=i (V'fcl· 
If we are given a frame T with frame operator G, then 

g = {G- 1 ' 2 \v) : \v) e ^} 

is a Parseval frame, see [TUl Theorem 4.2]. The projectors associated with Q 
form a POVM, since E^ G_1/2 l w ) H G ~ 1/2 = G^GG- 1 ' 2 = I holds. 

If we have a Parseval frame Q in C n with n 2 elements such that the asso- 
ciated projection operators are linearly independent and the frame elements 
satisfy the approximate symmetry (j3J), then the projectors corresponding to 
the frame Ç form an ASIC-POVM. 

2.3 Construction I: Pruning MUBs 

The first construction of ASIC-POVMs is based on the idea to select a suit- 
able collection of n 2 vectors from a set of n 2 + n vectors that form a maximal 
set of n + 1 mutually unbiased bases of C n . Our goal is to choose n 2 vectors 
such that the corresponding projection operators are linearly independent. 
We recali a known fact that belongs to the folklore of mutually unbiased 
bases; it is implicitly contained in [213 ES], and more explicitly in [22], and 
our proof is based on the latter. 

Lemma 3. Suppose that B a = {v a ^ : < b < n}, with < a < n, are n + 1 
mutually unbiased bases of C n , and let 

V = {Kb) (v*,b\ ■ 0<a<n,0<b<n} 

denote the associated set ofn 2 + n projectors. The n 2 projection operators in 
V* = {\v a ,b) ( v a,b\ '■ ( a , b) = (0, 0) or b ^ 0} C V are linearly independent. 
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Proof. First, suppose that a linear relation 



n n—1 



^2 7*,b Kb) (Va,b\ = (6) 
a=0 6=0 

holds for some G C. We are going to show that this has some rather 
strong consequences for the coeficients 7^. If we apply the projection oper- 
ators from V and take the trace, then we obtain a linear system of equations 
Ag = 0, where 

A = [trflüo.b) (ü aj6 | \v c4 ) (V c ,d\)}(a,b),(c,d) 

and g = (j aí b) is a column vector. The matrix A is block-circulant, 



.4 



( In n Jn " • • n Jn \ 

ij J ... Ij ij 



. I / 1/ ... 1/ 7 ; 



with nxn identity matrices in the diagonal blocks, and múltiples of the nxn 
all-one matrix in the off-diagonal blocks. 

If we subtract two equations in Ag = that belong to the same block 
row, then we find that 7^ = 7^ holds for all < b, d < n and all indices a. 
Therefore, the coefficients 7^ do not depend on the value of b. 

Finally, suppose that the left hand side of © consists of a linear com- 
bination of projectors belonging to the set V*, meaning that the coefficients 
7a,6 = when a ^ and 6 = 0. It follows that 7^ = holds whenever 0, 
since ^ a ,b = 7a,o by our previous observation. Therefore, the left hand side 
of (jïïj) reduces to 



n—1 n—1 



^2lo,b \v ,b) (v ,b\ = 5Z 7o '° \ v °^ ^°' fe l = °- 

6=0 6=0 

Thus, we must have 7 0i o = 0. Therefore, we can conclude that the projectors 
in P* are linearly independent, as claimed. □ 

We also recali the bàsic construction of MUBs in prime power dimension; 
see, for instance, |27| Hüj. 
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Lemma 4. Let q be a power of a prime p > 3. Define 

Wa,b) = q- 1/2 (ep (Tr Fq/Fp (ax 2 + bx))) x&q G C. 

Then the Standard basis Bo together with the bases B a = {\ip a ,b) '■ b G ¥ q }, 
a G F q , form a set of q + 1 mutually unbiased bases of <C q . 

Our first construction of ASIC-POVMs is given in the next theorem. 
Theorem 5. Let q be a power of a prime p > 3. Let 

\iP a , b ) = q- 1 ' 2 (e p (Tr ¥q/¥p (ax 2 + bx))) x&q G C q 

for all (a,b) G F, x and \ip afi ) = (6 a)X ) xe¥q for all a G F g . We define 
E a ,b = \^a,b) (i>a,b\ /q and 

G = YY1 Ea < b - 

Then the set {F a , b : a, b G F,}, with F a , b = G~ 1/2 E a , b G~ l/2 , is an ASIC- 
POVM. 

Proof. The linear independence of the operators E ajb follows from Lemma 
It remains to show that the matrix G _1 is close to the identity and that F ajb = 
G- 1/2 \ip a , b ) {ip a ,b\ G- 1/2 indeed forms an ASIC-POVM. First, we explicitly 
compute the frame operator G. We have 

G = - ^ \lp a ,b) (lpa,b\ +~Iq 

q q 

(a,fe)GF 9 xFj 

= 4 ( Y e p { Tr v q /F P (a( x2 -y 2 ) + K x ~y))) \ x ) (y\ 

1 \ x,y,a,b&¥ q 

- Y e p (^/wMi^-y 2 ))) \ x ) (v\ ) + 

x,y,ae¥ q ) 1 

We notice that (x | G\ x) = 1 for x G F q , (x\G\—x) = —1/q for x G F*, 
and (x\G\y) = for x, y G F q with x ^ ±y. Therefore, we can express the 
operator G in the form 

G = I a --Q + - 10) (01 , where Q = V Ix) (-x\ . 
q q " 

* H x& q 
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Using the structure of G, it follows that the inverse is given by 



G- 1 = 1 + 



i q + -^Q - |0) (0|. 



q 2 — 1 J q 2 — 1 q — 1 

Observe that the set of vectors {iV'o.b) : ( a ; &) £ x } is hivariant 
under Q, since Q iV'a.b) = iV'a.-ft)- Recali that by Lemma |U the bounds 

1(^,61^)1 <g" 1/2 

hold, whenever (a, b) ^ (c,d). Defining \tp a ,b) = G~ x l 2 \ip a ,b) we get 
F a , b = G~ 1/2 E aíb G~ 1/2 = l^iM/q 

and obtain that 

q 2 tr(F atb F c4 ) = 



< 



2 



1 



+ 



q 2 -l 

q 



q 2 -ï 



< - 1 + 



1 



(l/laAlpc-d) 



1 

+ - 



q 2 — 1 J q \q 2 — 1 



q(q - K 

2 



q(q - 1) 



1 + 0(1) 



This shows that the rank-one operators F a ^ form an ASIC-POVM. □ 



2.4 Construction II: Perturbing MUBs 

We now describe a second, different, method to obtain a set of n 2 vectors 
such that the corresponding projectors span the space of all n x n matrices. 
This construction of ASIC-POVMs works for all dimensions n such that n is 
an odd prime number. 

We note that all arithmetic operations in any expression involving ele- 
ments of F p and real numbers are performed over the real numbers (where 
each element of ¥ p is represented by an integer in the range [0, p — 1]). For 
example, for a real rel and a,x G F p , the power r ax means r u , where the 
integer u = ax can be of size (p — l) 2 . 
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Theorem 6. Let p be an odd príme number, and let < r < 1 be a real 
number. For a, b G F p define 

/ 1 — r 2a 

l^> 6 > = V i_ r 2 P a {r ax e p {ax 2 + fec))*^, G C p . 

and let E a ^ = \<p a ,b) (¥W>| /P- ^ en ^ e ^a,6 ° re linearly independent. Fur- 
thermore, let 

a,&GF p 

T/ien /or r = 1 - p" 3 the set {F a ^ b : a, b G F p }, wift F a , 6 = G~ 1/2 E a)b G~ 1/2 
is an ASIC-POVM. 

Proof. First, we show that the matrices E a ^ are linearly independent. Note 
that instead of considering the normalized vectors \<p a ,b) it is possible to 
consider the vectors \<f a ,b) — (r ax e p (ax 2 + bx)) xe ^ and to show that the 
corresponding projectors E a>b = \<p a ,b)(ípa,b\ are linearly independent. 

We use the technique introduced in [Tïï| to check whether the matrices 
E a fi are linearly independent. To each nxn matrix E ayb we associate a state 
\E a ^) which is simply the row-wise concatenation of the entries of E ajb as a 
vector of length n 2 . Then the matrices E ajb are linearly independent if and 
only if the matrix S = - J2 ab( z ¥p \E a ,b)(E a ,b\ has full rank. We obtain 

P S=J2 \Ea,b){É a , b \ 
a,be¥ p 

= r< x+ y +u+ ^e p {a{x 2 -y 2 -u 2 +v 2 )) 

a£F p 
x,y,u,v£ Fp 

x e p {b{x— y— u+v )) \x,y) (u,v\ 
beVp 

=p r a ^ +u+ ^e p (a(x 2 -y 2 -u 2 + v 2 ))6 x _ y ^ v \x,y}( Uj v\. 

aÇFp 
x,y ,u,vÇ: Fp 

The rows of S are labeled by pairs (x, y), with x, y G F p , and the columns 
by pairs (u,v), with u, v G F p . We first note that S can be written as a 
block-diagonal matrix of p submatrices, each of size p x p, if the rows and 
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columns of S are suitably rearranged. For < i < (p — 1), we define the 
sets Ci = {(x, x + i) : x G ¥ p } (we treat i as an element of ¥ p so a; + i is 
computed in ¥ p too). We now order the rows and columns according to the 
list C = IJf=o With respect to this basis we obtain that 



S = A © Ai © ... © A p _i, 



with pxp matrices 



Ai 



- x,u£¥ p 



where y = x + i and v — u + i. Hence, we obtain that 



^ r 2a{x+u+i) ep (2ai(u-x)) 



a£F p 



x,ue¥ p 



and have to show that this matrix is invertible for < i < (p — 1). In order 
to do so, we first observe that x \— > x — i/2 defines a permutation of the rows 
of any p x p matrix and that similarly u \— > u — i/2 defines a permutation of 
the columns (hereafter i/2 is computed in ¥ p ). Applying both the row and 
the column permutation to Ai we obtain the matrix 



Bi = 



^ r 2a(x+») ep (2ai(u-x)) 



a£¥ v 



x,u£¥ p 



Note further that x i— > x/2 and y ^ y/2 induce permutations of the rows 
and columns of any p x p matrix. Applying this to Bi we obtain the matrix 



r a(u+x) e p (ai(u - x)) 

a£¥ p 



- x,u£¥ p 



^(r^e p (-^)) a (r«e p (m)) a 

a&p 



15 



Since |r e p (z) 



\re p {-i) 



r and < r < 1, by the property of Vander- 
monde matrices [32] we conclude that C, is invertible for all i = 0, . . . ,p — 1 
which implies that the matrices Aj, and finally S are invertible. Arguing 
as in the proof of Theorem|Slwe have established the informational complete- 
ness of the projectors corresponding to the normalized vectors \<f a ,b)- 

Next, for r = 1 — p~ 3 we derive the bound 



\(<fa,b\fc,d)\ < 



1 + 0(1) 



(o,6) 7^ (c,d). 



(7) 



We have 



\(<Pa,b\ ( Pc,d)\ 



1—r 



2<i 



1—r 



2c 



1 r^2pa V ^ j~2pe 



J2 r 



(a+c)x 



e p (ax 2 +(3x) 



(8) 



where a = a — c and (3 = c — d. 

We frequently use that r* = 1 + 0{t/p 3 ) for any í = 0(p 3 ). In particular, 



1—r 



2a 



1—r 



2c 



J j,2pa V ^ ^2pc 



-1/2 



-1/2 



E 

, xGF„ 



,2ax 



E 



,2ca: 



Furthermore, 



A a+C > e p {ax 2 +(3x) 



< 



< 



< 



y2 r 2px 

x£W p 

:i+o(i)K 1 . 



(p(l+0(l/p))) 



-1 



^(l + 0(l/p))e p (ax 2 +/3x) 



e p (o;x 2 +/3x) 



+ 0(1)<VP + 0(1). 



Substituting the above bounds in (jHJ, we derive (JZJ). 

Next, we note have that E af) = \<p atb ) {ip a>b \ jp and G = Y. a ,be^ P Ea ^ We 
now compute the frame operator G and show that G^ 1 is close to the identity. 
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Similarly to the proof of Theorem |5] this implies that F a ^ = G 1 l 2 E a f ) G 
indeed forms an ASIC-POVM. We have 



G 



( 2 \ 

i _ r 2a 

E v(1 - r 2 Pa) ra(x+v) *M* 2 -y 2 ) + K*-v)) l*> (y\ 



V 

x,y,&p 



V 

/ a \ 



/ 



VClfcJTp 



/ 



1 r 2a 



Recalling that r = 1 — p 3 , from the Taylor expansion we obtain that 

l _ r 2a l + o(l) 



J j.2pa 



P 



Hence 



aSF„ 



1 — r 



2« 



„2a:c 



j,2pa 



P 



aeF„ 



,2ax 



l + o(l) 

P 



(p + 0(l)) = l + o(l). 



Finally, we deduce that 

p 2 tï(F a)b F c4 ) = \((Pa,b\&c,d)\ 2 = KwlG^Vc.d)! 2 
= (l + 0(l))|(^ 0l6 |^c,<i>| 



2 1+0(1 



P 

which implies that the rank one operators F a ^ form an ASIC-POVM. □ 



3 Relaxed MUBs and SIC-POVMs 



3.1 Motivation 

The constructions in the previous sections required the existence of n + 1 
mutually unbiased bases in C n . Currently, it is not known whether such 
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extremal sets of mutually unbiased bases exist when n is divisible by two 
distinct primes. We show that approximately mutually unbiased bases exists 
in all dimensions. Furthermore, we show that if we slightly relax the ASIC- 
POVM condition, then we can obtain in any dimension systems of vectors 
that approximate SIC-POVMs. For these constructions, we need some results 
about the distribution of primes and bounds on exponential sums, which we 
summarize below. 

These constructions work in any dimension n but attain its maximal 
strength for n of a certain arithmetic structure. In particular, if n — p — 1 
for a prime p, then we show the existence of AMUBs in C n . It is not known 
whether maximal sets of MUBs in these dimensions exist. 

Besides the general mathematical interest in the vector systems derived in 
this section, we expect that our approach and the new technique introduced 
in the following will lend themselves to some further applications in quantum 
information processing. 

3.2 Analytic Number Theory Background 

First of all we recali a remarkable result of [3] on gaps between consecutive 
primes. 

Lemma 7. For any sufficiently large x, any interval ofthe form [x — x ' 525 , x] 
contains a prime number. 

We now need some bounds of exponential sums with polynomials. 

Let p be a prime number and let ¥ p be a field of p elements. We always 
assume that ¥ p is represented by the elements {0, . . . ,p — 1}. 

The following statement is a variant of the celebrated Weil bound, see 
Example 12 of Appendix 5 of [H| as well as Theorem 3 of Chapter 6 in [3*ü] 
and Theorem 5.41 and comments to Chapter 5 of [SI])- 

Lemma 8. Let x be a nontrivial multiplicative character of ¥ p of order s 
Suppose that G(X) e F P [X] is not, up to a nonzero multiplicative constant, 
an sth power in ¥ p [X]. Then for any polynomial F(X) G F P [X] of degree d 
we have 

£ e p {F{u)) X {G{u)) 

u=l 



<(d + Z/-l)j9 1/2 , 
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where v is the number of distinct roots of G in the algebraic closure of¥ p . 

We now use Lemma|ïï]to estimate some mixed exponential sums with two 
exponential functions. 

Lemma 9. Let F(X) 6 F p [Jf] be of degree d > 2. Then, for any integer k, 



^2 e p (F(u))e n (ku) 



u=l 



0(p 2 / 3 ), i/d = 2, 
0(p 3/4 ), z/rf>3. 



Proof. Because each term in our sum is of absolute value | e p (F(u)) e n {ku) \ 
1, for every integer v > we have: 
p p 

e P (F(u)) e n {ku) = ^ e p {F{u + v)) e n {k{u + v)) + 0(v). 

u=l u=l 

Thus for every positive integer m we have 

p m—l p 



m e P (F(u)) e n (ku) = e p( F ( u + v )) e n{k{u + v)) + 0( 

u=l 

Therefore 



m 



v=0 u=l 



P 



^ e p( F ( U )) ^n(ku) 



u=l 



< — W + 0(m), 
m 



(9) 



where 



W 



p m—l 



<- E 

u=l 



J2 Yl e p( F ( u + v )) e M u + v )) 

u=l v=0 
p m—l 

^ e n (ku) ^2 e p( F ( u + v ï) e n(kv) 

m—l 

Y e p (F(u + v))e n (ki 



v=0 



By the Cauchy inequality we obtain 



w 2 < P J2 



m—l 



Y e p (F(u + v))e n (ki 



u=l 
m—l 



v=0 



p e n (k(v - to)) e p (F(u + v) - F(u + w)). 



v,w=0 



u=l 
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We now examine the polynomial F VjW (U) = F(U + v) — F(U + w). By the 
Taylor formula we have, 

F V , W (U) = F(U + v )-F(U + w) = J2 ^-j — U v - 

Clearly F^(v) = F^ (w) is possible only for v = w. For such m pairs 
of v and w we estimate the sum over u trivially as p. Otherwise we estimate 
these sums as (d — 2)p 1 / 2 by Lemma|Hl getting 



W 2 



O (mp 2 ) , if d = 2, 

O (mp 2 + m 2 p 3/2 ) , if d > 3. 



Thus by Q 

v 

^2 e p (F(u))e n (ku) 

u=l 

Taking m = |"p 2 / 3 ] we conclude the proof. 



O [mr ll2 p + m) , if d = 2, 

O (m-y^p + p 3 / 4 + m) , if d > 3. 



□ 



We also need a special case of the classical Weyl bound which we present 
in the following form, see Lemma 3.6 of [2^j or Lemma 2.4 of [IH] for a similar 
statement in full generality. 

Lemma 10. Let F(X) G F P [X] be of degree d > 2. Then for any fixed e > 
and any integer h < p, 



u=l 



O h 



1+e 



1 1 D 
- — I h 7^7 

h p hr 



l/2 d 



3.3 Arbitrary Dimensions 

We now describe a construction of a vector system, which satisfies ()3]) exactly 
and also gives a certain approximation to In fact we are able to get n d + l 
(rather than n + 1) orthogonal bases with this property, where d > 1 is any 
integer. 

Let Td be the set of polynomials of the form 

d+l 

f{x) = Y. a » xr 

v=2 
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with integer coeficients in the range < a v < n — 1, v = 2, . . . , d + 1. Thus 
j^Td = n d . Let p be the smallest prime with p > n. For each / G Td we 
consider the basis 

B / = >■■■> IV>/.n)} where = (e p (/(u)) e„ (m))" =1 . 

v ^ 

(10) 

Theorem 11. For any integer d > 1, í/ie Standard basis and the n d bases 
Bf,f& Td, given by ÍÍTÜj) are orthonormal and satisfy also 

w/iere f,g eT d U {0}, f ^ g, l<i,j <n. 

Proof. The orthonormality of each basis follows from the identity 

1 - 

U=l 

Clearly, if / £ T d and g = (or / = and # G T d ) then = ™~ 1/2 - 

Thus it remains to estimate 

1 n 

for /, g G T A , f ^ g and 1 < i, j < n. 
Because 

I e„ (/(u) - g(u) + (i - j»| = 1 
and by Lemma[7|we have 

1 - 

- Y] e P Ui u ) - 9{u)) e n ((i - j)u) 

u=l 

= -J2 e p (f( u ) - 9(u)) e n ((i - j)u) +0(\p-n\) 

u=l 

= - f> p (/(«) - 9{u)) e n ((i - j)u) + 0(n-° 475 ). 

u=l 
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Hence, 



1 



-0.475^ 



u=l 



Because f(X) — g(X) is a polynomial of degree at least 2, Lemma |ïï] yields 



^ e p (f(u) - g(u)) e n ((i - j» 



u=l 



0(p 2 / 3 ) = 0(n 2 / 3 ), d = 1, 
0(p3/ 4 ) = 0(n 3 / 4 ), d > 2, 



which concludes the proof. □ 
As before, let p be the smallest prime with p > n. We consider n 2 vectors 

B = {\1> f ) : f E where \^ f ) = -±= (e p (f(u))) n u=1 . (11) 



Theorem 12. Leí p 6e the smallest prime with p > n. Then, the vector 
system B of n 2 vectors given by (fTTJ) is normalized and also satisfies 

\(^ f )\<(2 + O(n- 1 / 10 ))n- 1 / 2 , 
where f,g 6f 2 , / ^ 9- 
Proof. Obviously, we have 

1(^/1^)1 = 1, /Gf, 
Put h = p — n. Then for /, g e JF 2 , / 7^ g, we have 



K^hMI = - 



1 

n 



1 

n 



XI e p(/( n + u ) -g( n + u )) 



u=l 



e p(f( u ) - 9(u)) 

u=l 

p 

Yl e p(f( u ) - 9(u)) 

u=l 

By Lemmas 151 and ITÜ1 we have 

\(^ g \^f)\ < 2p 1/2 n- 1 + O (^h^n- 1 

= 2p 1 ' 2 n- 1 + O {n- 1 (h^ +£ + h 1+e p- 1 ' 4 + /i^+y/ 4 )) . 



1 1 p 

h p h 3 



1/4" 
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By Lemma[7|we see that h = Ofjo ' 525 ), therefore 

for sufficiently small e and sufíiciently large p. Noting that 

pV2 = („ + o (n a525 )) 1/2 = n 1 / 2 (1 + (n- a475 )) 1/2 = n 1 / 2 + O (n 02625 ) 
we finish the proof. □ 

3.4 Special Dimensions 

Here we give some improvements of the constructions of Section 13.31 for the 
vàlues of n for which the smallest prime p with p = 1 (mod n) is sufficiently 
small. 

Let p be the smallest prime such that p = 1 (mod n). Let X n be the set 
of n characters of order n modulo p and U n the subgroup of residues of order 
n in W* . In particular jftJ. n = n. It is know that X n is a cyclic group, so for 
some character \Çi X n all other characters of X n are given be the powers x\ 
z = 1, . . . n. 

For / G F P [X] of degree at most d and the above character x 6 í„ we 
define 

B/ = {|^/,i>.···.lfe>} where |^> = ^ (^(/(tOMu)*) • (12) 
Let be the set of polynomials of the form 

d 

f{x) = Y, a » xr 

with integer coefhcients in the range < a v < n — 1, u = 2, . . . ,d + 1. 

Theorem 13. For any integer d > 1, the Standard basis and the n d bases 
Bf, f G Gd, given by (I12j) are orthonormal and satisfy also 

w/iere /,o G J^U {0}, f ^ g, 1 <i,j <n. 
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Proof. We have, 

(^> = ;Ex(«)' 

We also have 

Now, using LemmalHl we conclude the proof. □ 

Corollary 14. Let p be a príme and let n = p — 1. Then an AMUB exists 
in dimension n. 

Proof. We apply Theoremll3lfor p = n + 1 and for d — 1. Hence, we have 
n 2 + 1 orthonormal bases such that inner products between their components 
are all bounded by (n + l) 1 / 2 ?^ 1 = n~ x l 2 + 0{rT x ). □ 

4 Remarks and Open Qüestions 

The qüestions about finding SIC-POVMs and MUBs can be reformulated as 
a spherical design question in the vector space C n (see Zauner's thesis jSO] 
and also [2H1 HI] ) • Thus it is possible that the techniques of |2H1 , as well as of 
more recent works, see a very inspiring survey [SS], rnay apply to the problem 
of constructing systems of n 2 equiangular lines in C n , that is, SIC-POVMs. 
In fact, it is quite possible that with some adjustments they may also apply 
to MUBs. 

It is widely believed, see |H], but remains unproved, that the smallest 
prime p = 1 (mod n) satisfies the bound 

p = 0{n log 2 n). 

In this case, the bound of Theorem ITÏÏ1 becomes 0(n _1//2 logn). Thus, it is 
quite possible that the construction of Section EOl is always superior to those 
of Section 13.31 



r i, i = 

\ 0, 
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Finally, we remark that many of the results of this paper remain un- 
changed if one uses prime powers q = p r (and thus general finite fields ¥ q ) 
instead of just primes p. In particular, Corollary El holds true in this more 
general setting. Hence, in summary, we have shown that ASIC-POVMs and 
AMUBs exist for any prime power dimension q. Moreover, we have shown 
that approximate versions of mutually unbiased bases and SIC-POVMs exist 
in any dimension if we are slightly more liberal about our constraints on the 
angles. 
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